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Abstract 



A nilpotent quotient algorithm for finitely presented Lie rings over Z (LieNQ) is 
described. The paper studies the graded and non-graded cases separately. The 
algorithm computes the so-called nilpotent presentation for a finitely presented, 
nilpotent Lie ring. A nilpotent presentation consists of generators for the abelian 
group and the products expressed as linear combinations for pairs formed by gen- 
erators. Using that presentation the word problem is decidable in L. Provided 
that the Lie ring L is graded, it is possible to determine the canonical presen- 
tation for a lower central factor of L. LieNQ's Complexity is studied and it 
is shown that optimizing the presentation is NP-hard. Computational details 
are provided with examples, timing and some structure theorems obtained from 
computations. Implementation in C and GAP interface are available. 



1 Introduction 



The nilpotent quotient algorithm for finitely presented Lie rings — LieNQ as 
we shall refer to it in the sequel — operates with Lie rings over Z. Its essential 
goal is to compute a sufficient and useful presentation for a given Lie ring. It 
actually means to compute a presentation for its abelian group and to compute 
the structure constants in order to determine the Lie ring structure. By efficient 
and useful we mean that several important pieces of information can be read 
off immediately (e.g.: nilpotency, nilpotency class etc). Furthermore, by this 
presentation the so-called word problem is decidable, although it is known to 
be undecidable in the most general case, on word problem in groups and Lie 
algebras see (Hall, 1958). 

Such algorithms have existed for several decades for groups and Lie rings 
as well. Recall, that groups and Lie rings have similar structure if we think 
the commutator as a second operation in a group. It has a Jacobi-like property 
and an identity similar to the distributivity and it is anti-commutative, see for 
example (Robinson, 1982). This allows us to alter the known group algorithms 
for Lie rings. 

The first widely known (and maybe the most successful) algorithm is due 
to (Havas & Newman, 1980). It has several implementation in Fortran, C 
and GAP. The description of these were reported in (Celler, Newman, Nickel & 
Niemeyer, 1993) and (Havas & Newman, 1980). They were successfully used to 
compute the order of some Burnside groups. 

From our point of view, probably the most important algorithms are presen- 
ted in (Havas, Newman & Vaughan-Lee, 1990), (Nickel, 1993) and (Nickel, 
1995a). The algorithm described in this paper is a mixture of the latter men- 
tioned ones, however it differs at several points from those. 

There are three main sections in this paper. The first part contains the basic 
properties of the nilpotent presentation and describes an algorithm to compute it 
in the most general case. The next part is devoted to graded Lie rings. We shall 
see that in the graded case we can simplify LieNQ and, using other techniques, 
it is possible to determine the isomorphism type of the underlying abelian group 
of such a Lie ring. The last part contains computational issues. We shall see 
how the implementation works and some examples with timing will display the 
power of LieNQ in both graded and non-graded cases. 

The basic properties of Lie rings can be found in any textbook on that topic. 
We refer (Humphreys, 1972) for details. The first lemma is of fundamental 
importance for our goal, since it essentially claims that it is realistic to build a 
nilpotent quotient algorithm for finitely presented Lie rings. 

Lemma 1.1 In a finitely presented Lie ring each lower central factor is finitely 
presented. 

The proof can be done by a similar argument to that presented in (Robinson, 
1982) 5.2.6, but one might want to see (Schneider, 1996) for more details. 
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2 The Nilpotent Presentation of a Lie Ring 



Let L denote a Lie ring and L % the ith term of its lower central series. Product in 
L is denoted by brackets — as it is usual — and we use the left-normed convention 
that is, [a, b, c] = [[a, b],c]. If L is finitely presented Lemma 1.1 implies that 
L l I ' L l+1 is finitely presented for each i. Since at the present stage we are only 
interested in the nilpotent factor L/L k , wc think L is nilpotent of class k — 1. 
In other words we think L/L k is equal to L. In this case L has a so-called 
nilpotent presentation, that is a tuple (Q, S), where Q is a finite set of generators, 
{ai, . . . ,a r } say, and S is a finite set of relators of the form 

7i • a» = • ai+i H h a» ir ■ a r for some i, TR 

[aj, a»] = ■ aj + i + • • • + f}j t i tT ■ a r for each j > i, PR 

where the a, [3, 7 are integer coefficients, furthermore, 7^ > 1 where applicable. 
Note that the relations of type TR are sometimes referred to as torsion relations, 
while those of PR are known as product relations. A more detailed description 
and the proof for existence of such presentation can be found in (Schneider, 
1996), while it is easy to see, that a nilpotent presentation always defines a 
nilpotent Lie ring. Observe that we demand the existence of TR relations for 
each pair where j > i, but only some of those of type PR (even possibly 

none). Throughout this paper I denotes the set of indices, so that i £ I iff aj 
has TR relation. 

One easily sees that the relations of type TR determine the underlying abclian 
group structure for L, while the relations of type PR define its ring structure. 
The (3 can be viewed as structure constants for the Z-module defined by Q and 
the TR relations. 

For our purpose this concept will be too general, so we shall keep some 
restrictions. Define a weight function u : Q — > N as follows. Let u be an 
increasing function in i, such that the following holds. 

1. w(ai) = 1, 

2. if the nilpotent presentation contains a relation of the form 

[a ,ai] = wji, 

where Wji is a sum of integer multiples of generators from Q, then all 
generators a/j £ Q such that w(dfc) < uj(aj) + uo{ai) have coefficients zero 
in Wji. 

An other restriction, if G Q and u>(dh) > 1, we require au to have a 
definition, i.e. a relation of the form 

(1) [aj,ai] = a k 

where u>(ah) = w{a,j) + u>{ai) and u)(di) = 1. It might easily happen that a 
nilpotent presentation contains more then one relations of the form (1) for some 
i, in this case we arbitrarily choose one for the definition for dj. 
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From now on by a nilpotent presentation we shall always mean weighted 
nilpotent presentation which satisfies the condition stated above for the existence 
of definitions. 

Suppose that we are given a nilpotent presentation 

(2) (Q | 7i • Oi = Wi, i e I, [a,j,a,i] = Wji, 1 < i < j < r ). 

Let L denote the Lie ring presented by (2). Clearly (2) does not care about Jac- 
obi identities, so L is not automatically coincides with the freest non-associative 
anty-commutative ring holding the relations of (2) in it. However we can arrange 
a collection process in L. Suppose, we are given a word-sum I expressed in term 
of the elements from Q . The collection process consists of two main steps. In 
the first step we substitute the products by the right hand side of the appropriate 
PR relation. In the second step we reduce those coefficients that are not less 
than the corresponding 7$ in (2). 

After proceeding so, we get a normal form for I subject to (2). That is, 

(3) I = a\ ■ on H h a r ■ a r 

where < on < 7$ whenever i £ I. An arbitrary I might have more than one 
normal form of type (3). However uniqueness of (3) is a crucial point, if we want 
to attack the word problem in L. We might want the nilpotent presentation to 
possess this property according to the following definition. 

Definition 2.1 A nilpotent presentation is said to be consistent if every element 
of L uniquely has a normal form. 

In Section 3.3 we shall see how to check if an arbitrary nilpotent presentation 
is consistent and develop a method to make it consistent if it is not. In general, 
it is an equivalent condition to that of Definition 2.1 that the element uniquely 
has a normal form. 

From the computational point of view it is important to see that all the 
relations of type PR are not necessary to determine the structure of L. Our first 
result is devoted to this observation. 

Lemma 2.1 Let L be a Lie ring given by the nilpotent presentation 

(Q \li ■ ai = Wi, i E L, [a j ,a l ] = w jl , 1 < i < j < r ). 
Then L is already determined by the presentation 

(4) ( Q I 7i • a l = Wi, i e L, [aj,ai] = w oi , 1 < i < j < r, uj(a l ) = 1 ). 

Proof. Essentially we want to prove that all PR relations can be expressed by 
those listed in (4). We proceed by an induction argument on w(oi). If w(aj) = 1 
then we arc done. Suppose now, that w(aj) = k > 1. Using the definition for aj 
and the identities of a Lie ring one has: [aj,aj] = [aj, [a,k,ai]} = — [a^, a;, a.,] = 
[a;, aj, ak] + [aj, a^, ai] — [aj, a^, ai] — [aj, aj, a^]. And the latter ones are known 
by the hypothesis. In the first equality we substituted ai by its definition. □ 
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3 Computing the Nilpotent Presentation 



The aim of this section is to show how it is possible to compute a consistent 
nilpotent presentation for a nilpotent factor of a a given finitely presented Lie 
ring L and determine an epimorphism from L onto that nilpotent factor. 

The basic ideas — besides being simple — are similar to those used in 
(Havas, Newman & Vaughan-Lee, 1990) and (Nickel, 1995a). In fact, it is based 
on an induction argument. Suppose, we are given a consistent nilpotent pres- 
entation for L/L l and an epimorphism from L onto L/L l , we shall show how to 
extend them to a consistent presentation and an epimorphism for L/L l+1 . 

As computational tools, matrices over Z play an important role. Computing 
the row Hermite normal form is a crucial point of the algorithm. Since its 
details can be found in (Nickel, 1995a) (Schneider, 1996), (Sims, 1994), we do 
not emphasize them further in this paper. 

3.1 Computing the Abelian Factor 

Suppose that we are given a Lie ring L with finite presentation ( X \ TZ ) , where 
X is a finite set of generators, while TZ is a finite set of relators. For simplicity, 
suppose X — {xi, . . . , x n }. Set Q = {a\, . . . , a n }. Construct a map 4> : X — > Q, 
by simply saying 

Xi4> = gi- 

Furthermore we set all PR relations to be trivial and 7 = 0. We call the pres- 
entation constructed above a trivial nilpotent presentation and it is clear that it 
represents the free abelian Lie ring on n generators. It is straightforward to 
see that it contains LjL 2 as a factor ring. Indeed, let X be the ideal generated 
by the set TZ<j>, where we evaluate <j) as if it was a homomorphism elementwise 
on TZ. Then we have 

Lemma 3.1 LjL 2 is isomorphic to JF° b /Z. 

Proof. Both Lie rings are generated by n elements and satisfy exactly the same 
relations. It checks the isomorphic property □ 

What we do in practice is the following. Evaluate the relations in as 
described above and put the images in a matrix M. Then compute the row 
Hermite normal form M H for M . It is well known that the rows of M generate 
the same subgroup, in the free abelian group, as those of M H . If the first non- 
zero element, the zth say, of a row of M H is 1, i.e. we have a row of the form 

(0, . . . ,0, l,m i+ i, . . . ,m„) 

we simply throw Oj out of Q and modify the map <fi 

Xi4> = -m i+1 ■ a l+1 - ■ ■ ■ — m n ■ a n . 

Otherwise, if the leading element is greater then one, i.e. we have a row of the 
form 

(0, . . . ,0,to; > l,m i+ i, . . . ,m n ), 
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we introduce a TR relation 



mi ■ a.i = — mj+i • a l+ i m n ■ a n . 

and put i into /. 

Proceeding each row as described above, we obtain a presentation for L/L 2 . 
We set the surviving generators of weight one and it is straightforward to see 
that the presentation obtained so far is indeed a presentation for the abelian 
factor and <f> extends to an epimorphism. We summarize this result in 

Proposition 3.1 The so obtained presentation is a consistent nilpotent pres- 
entation for L/L 2 and <fi extends to a epimorphism from L onto L/L 2 . 

The untouched images under epimorphism will be referred to as definitions 
of generators of weight one. This terminology will be important later on. 

3.2 Extending the Presentation 

Now we shall see how the induction argument works in details. Suppose, we have 
a nilpotent presentation for the factor ring L/L 1 , an epimorphism from L onto 
L/L 1 . Extending the presentation will consist of three steps, i.e. 

1. extending the epimorphism </>, 

2. modifying the TR relations, 

3. modifying the PR relations. 

Extending the epimorphism means introducing a new generator for all Xi G X, 
where Xi<f> is not a definition. In other words, if Xi4> = Wi, we modify Xi<p so that 
Xi<j> = Wi + U. 

The torsion relations will change so that if i 6 /, 7, • a, = Wu, then we alter 
7i • Oi = wu + tu. 

Those of the PR relations that are not definitions are modified in a similar 
way. If [a,j, at] — Wji j > i and [dj, a*] is not a definition, let the new relation be 

[dj,d l ] = Wji + tji. 

Note, that all newly introduced generators U,tu,tji are different from one 
another. We introduce PR relations so that all of them are central. At this stage 
we do not alter /. We prove the following 

Proposition 3.2 The so extended presentation contains L/L l+1 as a factor 
ring. 

Proof. Since we did not alter the definitions, every relation holding in the exten- 
ded presentation is satisfied in L/L l+1 . □ 

In practice we use Lemma 2.1. This lemma essentially says, that we do not 
need to introduce new generators for all PR relations, but only for those of the 
form [dj,di] = Wji for j > i and dj is of weight one. The remaining ones can 
be computed as we have seen in the proof of this lemma. This computation is 
usually referred to as computing the tails. 
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3.3 Enforcing Consistency 



Recall, we called a nilpotent presentation consistent if each element in L uniquely 
had a normal form. The extended presentation investigated in the previous 
subsection might happen to be not consistent, as the case is in general. 

In order to investigate the uniqueness property, we introduce an operation 
on the nilpotent quotient. First of all let A be an abelian group generated by 
Q and the TR relations of the extended presentation. Introduce a map such 
that, 

{Wji if j > i, 
-Wji if j < i, , 
ifj=i 

provided that the nilpotent presentation possesses a PR relation of the form 
[aj, a,] = Wji for j > i. Then the following is true. 

Proposition 3.3 ip can be extended to a binary operation on A iff 

n 

Oj-[aj,ai}= ^2 <Xjk[ak,ai\, CI 
k=j+i 

where j £ I and the relation corresponding to dj is of the form 

Oj ■ aj = djj+idj+i + • • • + a m ■ a n . 

Proof. It is immediate that the condition is necessary. To see that it is sufficient 
as well, observe that the two expressions, w and w' say, are equal, it means that 
they can be transformed to each other by using TR relations. The condition 
of the proposition claims that such a transition respects the equality ip(w,v) = 
ip(w',v) for any v. Similar result is obtained in the second variable of ip by 
swapping the two variables. □ 

If A is equipped by that operation, it happens to be a non-associative ring. 
We prove the following 

Lemma 3.2 The nilpotent presentation of L is consistent iff A is a Lie ring. 

Proof. Observe, L = A/ J where J is the ideal generated by all instances of 
Jacobi identity. A normal word I represents the zero element in L iff I £ J. If 
A is a Lie ring then J = 0, so I must be the trivial word. On the other hand, 
if the nilpotent presentation is consistent, then all instance's of Jacobi identity 
collect to the trivial word, so J = and A = L holds. □ 

One can easily check the Jacobi identity, since it is a multi-linear identity, 
thus it is enough to check it on the triples formed by the abelian group generators. 
However we can restrict ourselves even more, by the following 
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Lemma 3.3 A nilpotent presentation (G, S) of a Lie ring L is consistent iff the 
Jacobi identity holds for each triple (aj, aj, a k ) G Q x Q x Q where 1 < i < j < 
k <n and <Xj is of weight one. In other words 

[di, a 3 ,a k ] + [a J ,a k ,a l ] + [a k , a i: aj] = C2 

for l<i<j<k<n and aj is of weight one. 

The above lemma appeared in its first form in (Vaughan-Lee, 1984) for pc- 
presentation of p-groups. It was modified for Lie algebras in (Havas, Newman 
& Vaughan-Lee, 1990), where one finds its proof. The Lie ring case can be 
handled in the same way as Lie algebras. 

3.4 Enforcing the Denning Relations 

One can easily observe that the presentation obtained so far is a presentation 
of the freest Lie ring of class I that happens to contain L/L l+1 as a factor ring 
in it. What is still needed is to assure that the images of the elements of 1Z, 
under the epimorphism, are 0. It simply means taking the factor ring over the 
ideal generated by the epimorphic images. By the induction hypothesis those 
images vanish in L/L l , so they must lie in L l /L l+l and we are allowed to use 
abelian group methods again, such as Hcrmitc normal form and other integer 
based strategies. 

In practice we proceed as follows. Put all elements obtained from the con- 
sistency relations (CI), (C2) and the defining relations together in a matrix 
M. Compute the row Hermite normal form M H for M. Using those rows of 
M H whose leading element is equal to 1, we can eliminate some of the generat- 
ors from the nilpotent presentation. The other rows express linear dependence 
among the generators over Z. They can be viewed as TR relations and added to 
the presentation. Extend the weight function to the newly introduced generators, 
by saying that they are of weight i and add the the indices of generators that 
were provided with TR relations to /. The following lemma makes sure that 
generators of weight i have definitions. 

Lemma 3.4 If L is a Lie ring and suppose that L/L 2 is additively generated 
as 

L/L 2 = (£ 1+ L 2 ,...,£ k + L 2 ) 
and L 1 ^ 1 /L l is additively generated as 

L i - 1 /L i = (£ 1 ' + L i ,...,£i' + L i ) 
then L l /L t+1 is additively generated as 

L i/ L i+i = ([£' m J n ] + L l+1 \l<m<l,l<n<k). 

The proof of the above lemma can be done essentially in the same way as in the 
group case. That proof can be found in (Sims, 1994) Proposition 2.6. But see 
also (Schneider, 1996) for a detailed study of the Lie ring case. 
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4 The Graded Case 



4.1 Simplifying LieNQ in Graded Lie Rings 

The following section is devoted to graded Lie rings. Gradedness of a Lie ring 
is denned in accordance with the following definition. 

Definition 4.1 Let L be a Lie ring, L is said to be graded if L splits into a 
direct sum of its lower central factors. 

Since a Lie ring of that kind possesses a relatively easy structure, we are 
allowed to simplify our algorithm in this case. 

First of all, we need not extend the TR relations any more since if dj G 
L l /L l+1 for some i £ I implies 7, • ai G L l /L l+1 , so the right hand side of the 
relation corresponding to a, cannot contain generators of weight greater than I. 

Similar thing happens when we extend the PR relations. At the Ith step we 
only alter the relations of the form 



We regard these facts when we introduce new generators and compute the tails. 

Now it should be clear that enforcing the consistency simply means enforcing 
the consistency relations (C2) for the triples (a,,Oj,afc) G Q x Q x Q, where 
l<i<j<k<n u)(ai)+u)(a,j)+u>(ak) = I and w(a,) = 1. And the relations of 
type (CI) are checked for only pairs (ai,aj) G Q x Q, where ui(a,i) +tu(aj) = I. 
Note that I still denotes the nilpotency class of the current factor. 

We are allowed to simplify it further, namely we need to enforce only those 
relations of 1Z which have weight I. 

The natural question arises. How can one recognize a graded Lie ring by 
looking at its finite presentation? We only mention here the well known fact 
that a Lie ring defined by homogeneous relations is always graded. 

4.2 Another approach — Canonical Nilpotent Presentation 

In a graded Lie ring we shall construct the so-called canonical nilpotent presen- 
tation instead of nilpotent presentation defined in Section 2. The concept of the 
canonical nilpotent presentation fits the concept of canonical presentation for fi- 
nitely generated abelian groups. In the sequel we shall use additive notation for 
abelian groups and + to denote the binary operation in such structures. Recall 
that C n denotes the cyclic group of order n. 

Theorem 4.1 Let (A, +) finitely generated abelian group. Then A can be writ- 
ten as a direct sum of cyclic groups, i.e. 



where ki\k 2 ,. ■ .,k r -\\k r . The (5) form is called the canonical decomposition for 
A. The canonical decomposition is unique for an arbitrary such A. 



[a,j,ai] = Wji 



for i > j i + j = I. 



(5) 



a = c kl e • • • e c kr e © • • • © c< 
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The proof of the theorem can be found in (Suzuki, 1982) Theorem 5.2. 

In the following we define the canonical nilpotent presentation for a graded 
Lie ring. 

Definition 4.2 A Lie ring presentation is said to be canonical nilpotent pres- 
entation, if it is of the form 

(6) (a 1 ,...a n \c i -a l for i E I, [dj, a,] = ti) tJ forj>i), 

where the Wij are normal words as in (3). 

Remark 4.1 Using the results of Section 3 we shall see that each finitely presen- 
ted graded nilpotent Lie ring possesses such a presentation. Essentially it is 
enough to see, that the generators corresponding to a certain lower central factor 
can be chosen so that the torsion relations are trivial in (6). Conversely, it is 
rather easy to see that a Lie ring defined by a canonical nilpotent presentation 
is always a finitely presented, nilpotent and graded. 

Remark 4.2 Actually we shall construct more than a simple presentation defi- 
ned in (6). We shall construct a presentation in which the lower central factors 
are canonically presented as abelian groups, according to Theorem 4-1- That is, 
we shall enforce the appropriate divisibility conditions for the Ci . 

The following theorem is an easy consequence of Theorem 4.1 and Definition 

4.2 

Theorem 4.2 Requiring the divisibility condition for each lower central factor 
in (6), the coefficients ji are uniquely determined by the isomorphism type of 
L. 

4.3 Integer Matrices 

In this section we collect the background knowledge needed to construct our 
nilpotent quotient algorithm. First of all, a matrix is called integer matrix if its 
entries are rational integers. 

Definition 4.3 Let A and B m x n be integer matrices. A and B are called 
equivalent if there exist P m x m and Q n x n integer unimodular matrices so 
that PAQ = B. 

Remark 4.3 The unimodular matrix P corresponds to elementary row oper- 
ations, while Q corresponds to elementary column operations. An elementary 
row (column) operation is negating a row (column), adding an integer multiple 
of a row (column) to an other row (column), swapping two rows (columns). 

Remark 4.4 A trivial but rather useful observation is that the Smith normal 
form of the matrix (ai, . . . a n ) is (gcd(ai, . . . , a„), . . . , 0). 

The following result can be traced back to (Smith, 1861). 
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Theorem 4.3 For each M mxn integer matrix there exists an S mxn diagonal 
matrix, such that the non-zero entries of S are non-negatives and Si_i,i_i 
for 2 < i < min(ra, n). S is called the Smith normal form for M. 

Remark 4.5 Theorem 4-1 and Theorem 4-3 look very similar. Indeed, they are 
closely related. It is not hard to see if we put the defining relators of an abelian 
group into a matrix, then abelian groups with equivalent relation matrices are 
isomorphic. Theorem 4-.1 has a standard proof as follows. Put the defining 
relations of A into a matrix M. Then compute the unique Smith normal form 

5 for M and read off the canonical presentation for A from S. 

(Hartley & Hawkes, 1970) provides a naive algorithm to compute Smith 
normal form for a given integer matrix. However, in practice it does not seem to 
be widely applicable. Besides the fact, that the intermediate entries are growing 
very fast, we shall need to compute the unimodular matrix Q too. Q is not 
unique and we might want to choose a capable Q with moderate entries. The 
above mentioned algorithm does not fit this goal as examples will point it out 
later. 

Werner Nickel (Nickel, 1995b) suggested the following algorithm. 

Algorithm 1 

Input: mxn matrix M 
Output: Smith normal form for M 
while M is not diagonal 
begin 

M := Hermite normal form of M 

M := Transpose of M 
endwhile 
Enforce divisibility on M 's diagonal 
return M 
end 

Lemma 4.1 Algorithm 1 computes the Smith normal form in finitely many 
steps. 

Proof. If after a transposition the first entry of the ith column divides all 
entries in the its column, we can reduce that column without producing trash in 
the upper half. 

If it does not divide some entry in its column, after a reduction the first 
element will become a smaller positive integer. After finitely many steps it will 
divide all column entries (in the worst case it becomes 1). □ 

More sophisticated algorithms exist to compute Smith normal form as repor- 
ted in (Havas & Majewski, 1994), (Havas & Majewski, 1995) and (Havas, Holt 

6 Rees, 1993). They use various techniques such as LLL, modular methods and 
pivoting strategy. They behave very well in practice, however, our matrices are 
rather sparse and Algorithm 1 is suitable for our purposes. 
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Example 4.1 In (Sims, 1994) °n Page 379 one finds an example with detailed 
analysis of its Smith normal form computation. Using the naive algorithm it 
was hopeless to compute Smith normal form for the matrix. Using Algorithm 1 
one can compute Smith normal form, producing 5276851227 as the maximum 
entry of Q. 

Algorithms with more sophisticated background do not behave much better. 
Its Hermite normal form was computed using ideas described in (Havas, Holt & 
Rees, 1993), implemented by B.Majewski in GAP. Norm driven pivoting strategy 
produced 3452231215618892621983848708913689 - 3.4- 10 34 as maximum entry 
in the transformer matrix. Their LLL based algorithm computed the Smith 
normal form with 298901 as maximum entry of P. 

Example 4.2 The other example arose from the following Lie ring presenta- 
tion. 

(x,y\ [y,x,x,y], [y,x,x,x] + [y,x,y,y]) 

Computing its 10th factor, LieNQ produced an 30 x 16 matrix. Naive approach 
implemented in C, with 32 bit arithmetic, produced integer overflow. Algorithm 1 
provided a transformer matrix with largest entry 12. Norm driven Hermite nor- 
mal form computation produced 62, LLL based algorithm produced 5 as largest 
entry in P. 

The idea behind Algorithm 1 is that we keep the number of column operations 
as low as possible. It seems to be a good idea using LLL based algorithm, 
however, it has worse complexity than that of traditional strategies. 

As we saw the main problem is to keep magnitude of the entries of the 
intermediate matrices and the transforming matrix under control. To find the 
optimal strategy — even in the simplest case — is hopeless, as stated in (Havas & 
Majewski, 1994). 

Proposition 4.1 Let us given a — (a 1; . . . , a n ), a vector consisting of n positive 
integers. To find a vector x, such that x is shortest either in Lq or in L^ norm 
and solves the equality x ■ a = gcd(al, . . . , a n ) is HP-hard. 

The following proposition is an easy modification of Theorem 3 in (Havas & 
Majewski, 1994). 

Proposition 4.2 Let us be given an integer vector a = (ai, . . . , a n ). The task of 
minimizing the vectors in the unimodular matrix Q, where 
aQ = (gcd(ai, . . . , a„), . . . , 0), with respect to the maximum norm (MINGCD), 
is NP-hard. 

For the definition and basic properties of NP-hard problems we refer (Davis 
& Weyuker, 1983). 

Let the abelian group A = (xi, . . . x n \ r\, . . . , r m } be given by generators 
and relations. As it was mentioned above, the canonical presentation can be 
constructed for A by computing the Smith normal form for the relation matrix 
of A. We emphasize this idea in details now. 
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The relation matrix is constructed as follows. If n = a^ixi + • • • + a,i^ n x n let 
M = (dij) the m x n matrix. Suppose that the Smith normal form for M is 







\ 



s = 



V o 



0/ 



Then 



(7) A — Csn © C S 22 © ' ' ' © C Srr (B Coo © ' ' ' © Coo 7 

where the number of infinite factors is n — r. Since S — PMQ for suitable 
unimodular matrices P and Q, S can be obtained from M by elementary row 
and column operations. Since a row operation does not alter the generating set of 
A, P can be forgotten, while Q can be used to obtain the isomorphism between 
the two presentations of A. 

Proposition 4.3 Let A, M, Q and S be as above. Suppose moreover that the 
cyclic factors in (7) are generated by a\, 02,- ■ -,a n respectively. Then 



( 01 ^ 




( ^ \ 


and 


Q 


( a\ 


\ a„ J 




\ x n ) 






\ a n 



Proof. See (Sims, 1994) Proposition 8.3.1. □ 



4.4 Computing the Canonical Nilpotent Presentation 

After this preparation we are ready to describe, how one can compute canonical 
nilpotent presentation for a finitely presented nilpotent graded Lie ring. 

Recall, we computed the matrix consisting of consistency relations and de- 
fining relations. We computed its Hcrmitc normal form. Now we shall compute 
its Smith normal form together with the transformer matrix Q and its inverse 
Q^ 1 . Using Proposition 4.3 one can compute new generators for the canonical 
presentation. Put, Q = (qij) and Q^ 1 = (qij). 

The condition for the existence of the definition cannot be kept any more in 
that case. We use the following approach instead. 

Suppose, at a certain layer we introduced new generators x\, ... ,x s , so that 
they correspond to the following products: 

Xi = [bii, &2i] where u)(bu) = 1 and uj(b 2 i) — u)(xi) — 1, 1 < i < s. 
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Using Proposition 4.3 one gets the generators oi,. . .,a s of the canonical presen- 
tation for that layer. Moreover one has that 







( [611,621] \ 




Hi 








V [6i s ,6 2s ] / 



So, we can introduce definitions for the canonical generators of the following 
form 

ai := &,i [611,621] H 1- 9i, s [6i 5 ,6 2s ] for 1 < i < s. 

That is, the definitions are linear combinations consisting of Lie products with 
generators of smaller weight. 

In Section 2 we already computed products of current weight. They need 
one more look as follows. After tail computation, we expressed those products 
in terms of newly introduced generators x\,...,x s . Now, those generators are to 
be replaced by the canonical generators 01,. . .,a s . To do that we use Proposition 
4.3 again. That is, 

Xi = qi,ia,i H h q lyS a s for 1 < i < s. 

We proceed the old generators occuring in the cpimorphic images in the same 
way. 

Lemma 2.1 is still true in the canonical case and its proof can be done in 
essentially the same way as presented before. Lemma 3.3 seemingly has no 
connection to the form of the definitions. However, having a deeper look at its 
proof, one finds that it heavily depends on the concept of definitions, but can be 
modified to the new approach too. 

The canonical nilpotent presentation computed so far depends on the trans- 
former unimodular matrix Q. Among many possible matrices, we want to choose 
the one with the smallest possible entries, because they will appear in the pres- 
entation. This is a hard problem as stated in 

Theorem 4.4 Let L be a finitely presented nilpotent, graded Lie ring. The task 
of minimizing the vectors in the right hand side of the epimorphic images and 
the PR relations, according to the maximum norm, (MINLIE) is NP-hard. 

Proof . We reduce MINGCD to MINLIE in polynomial time. Let us be given an 
instance (ci, . . . , c„) of MINGCD. Write down the finite presentation 

(xi, ... ,x n \ c\X\ H \-c n x n = 0,[xi,Xj]=0 for l<i<j<n). 

Clearly, this can be done in time proportional to n 2 . Use LieNQ on this pres- 
entation. LieNQ terminates at the abelian factor and produces the following 
output: 

Xi(fi = qi,\a\ H qi, n a n for 1 < i < n, 

(8) gcd(ci,...,c„)ai =0, 

{ai, a,j] — for 1 < j < i < n. 
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Where the transformer matrix Q = (qij). The ability of minimizing the vec- 
tors in (8) implies the ability if minimizing the vectors in Q and hence solving 
MINGCD. □ 



An other disadvantage of this approach is presented by the following theorem. 

Theorem 4.5 The absolute values of the coefficients appearing in the canonical 
nilpotent presentation cannot be bounded by a function depending only on the 
number of generators. 

Proof. Easy consequence of (Havas & Majewski, 1994) Lemma 4. □ 



5 Performance 



5.1 Implementing LieNQ in C 

The implementation for the LieNQ algorithm has been written in C program- 
ming language. The currently available version is the Version 2.0, which provides 
time measure, computations for graded and non-graded cases separately. 

The program takes a finite presentation (X,1Z) for a Lie ring L and the 
nilpotency class of required factor as its input. It outputs the TR, PR relations 
and the images of the elements of X under the cpimorphism. The nilpotency 
class can be omitted, in this case the program runs until the lower central series 
stabilizes or the algorithm goes beyond the capacity of the computer. 

Throughout the computation the data of the nilpotent presentation is stored 
in normal word form. At previous stages they were stored as coefficient vectors, 
but it turned out to be more efficient in an array consisting of structures, where 
the first element of a structure contains the number of a generator and the second 
one its coefficient (the idea is due to (Nickel, 1995b)). For instance, the Lie ring 
element 

• a\ + • a2 + 2 • 03 — 4 • 04 + • 05 

is stored as 

((3, 2), (4, -4), (0,0)). 

The first component of the last pair indicates an EOW (end of word) sign, while 
the second one, in fact, is arbitrary. 

The information of the nilpotent presentation is stored in several arrays as 
follows: 



Coef f icients [i] ji if i S I, otherwise, 

Power [i] Wi, right hand side of the TR relation 7^ • a^, 

Product [j] [i] wji, right hand side of the PR relation [aj, ctj] 
Epimorphism[i] the ith epimorphic image, 

Weight [i] u(a-i), the weight of a^, 

Dimension [i] number of generators of weight i, 

Def inition[i] the definition of a^. 
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Note that Definition [i] is an array of structures. A component of the 
array either contains the number of two generators and a positive integer, if the 
definition of dj is like a* = ^ °i ' [ a fc; , azj , or contains the number of the generator 
from X in the first component and in the second and the third if a, is defined 
as an epimorphic image. In the latter case the length of the array is, of course, 
one. 

The defining relations are stored as expression trees. See (Nickel, 1995a) for 
its description which emphasizes its advantages in details. 

During the computation graded and non-graded cases are distinguished ac- 
cording to whether the user switched on the -g option or not. There are differ- 
ent tails routines, that is Tails and GradedTails, and consistency routines, like 
Consistency and GradedConsistency, built in the code. Of course, introducing 
new generators is also a different task. 

It is, however, important to introduce the new generators in a right order. 
Lemma {commlemma makes possible to express all tails in terms of generators 
cifc of weight I where a,k — [aj , a,] for some a,j of weight I — 1 and a, of weight 
1. In order to use this fact in practice, we introduce those generators the last 
and when we compute the row Hcrmitc normal form for the matrix they will be, 
in fact, the only surviving generators and all the others disappear. In this way, 
we ensure that the new generators indeed have definitions, so the restrictions we 
kept are reasonable. In the graded case such restrictions are not kept. 

The main routines for matrix computation and the basics of the parser pro- 
gram are due to W. Nickel. The first mentioned one uses the GNU MP package 
to deal with long integers. 

LieNQ was successfully compiled on several Unix platforms, such that Linux, 
Free-BSD, SunOS, Solaris. 

5.2 GAP Interface 

LieNQ is might be invoked from GAP (Schonert et at, 1994) using a GAP 
interface. Since GAP does not know about Lie algebras until version 3.5, it is 
necessary to have at least GAP-3.5 to use the interface. We remark, that even 
GAP-3.5 knows only Lie algebras but not Lie rings over a ring, so LieNQ kills 
torsions in that case and GAP considers the result as Lie algebra over Q. 

The interface consists of three functions, two of them are auxiliary ones and 
one is worth to more interest. An auxiliary function is implemented to print a 
Lie ring presentation according to LieNQ's syntax, and one is implemented to 
print a usage-message in the case of possible errors. 

Now we describe how to use the function LieNilpotentQuotient. 

LieNilpotentQuotient(P, c) 
LieNilpotentQuotient(P) 

In the first case the function computes the canonical nilpotent presentation for 
the cth nilpotent factor of P. In the second case it computes the largest nilpotent 
factor of P if it exists. 
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In both cases P has to a be record consisting of at least two fields: generators 
and relators. A relator is arbitrary Lie element, using the lcft-normcd convention. 
In complex cases it is a good idea to give relators as strings, since addition is 
not defined between two lists. 

The function returns a table Lie algebra record containing a field for the finite 
presentation P, one for the epimorphism and one for the lower central series of 
the Lie algebra defined by the output presentation. 



6 Some Sample Computations 

In this section we present some examples, that was tried by LieNQ. We consider 
the following Lie ring presentations. 

L\ ={x,y) 

L 2 =(x,y\ [y,x,y], [y,x,x,x,x,x]) 

L 3 =(x,y \ \[[y,x],x,y], [[y,x],y,x], [[y,x],x,x] + [[y,x],y,y]) 

L A ={a,b,c,d,e | [b, a], [c, a], [e, c], [e,d], [d,a] = [c, b], [d,b] = [e, a], [d,c] = [e, b]) 

L 5 = (ei, e 2 , e 3 , e 4 , e 5 , e 6 , e 7 , e 8 , e 9 , e 10 

[e2,ei,ei], [ei,e 2 ,e 2 ], [e 3 ,ei], [e 4 ,ei], [e 5 ,ei], [e 6 ,ei], [e 7 ,ei], [e 8 ,ei], 

[e 9 ,ei], [ei ,ei], [e 3 ,e 2 ,e 2 ], [e 2 ,e 3 ,e 3 ], [e 4 ,e 2 ], [e 5 ,e 2 ], [e 6 ,e 2 ], [e 7 ,e 2 ], 

[eg, e 2 ], [e 9 , e 2 ], [ero,e 2 ], [e 4 , e 3 , e 3 ], [e 3 , e 4 , e 4 ], [e 5 , e 3 ], [e 6 , e 3 ], [e 7 , e 3 ], 

[eg , e 3 ] , [eg , e 3 ] , [eio , e 3 ] , [e 5 , e 4 , e 4 ] , [e 4 , e 5 , e 5 ] , [e 6 , e 4 ] , [e 7 , e 4 ] , [e 8 , e 4 ] , 

[e 9 ,e 4 ], [eio,e 4 ], [e 6 ,e 5 ,e 5 ], [e 5 ,e 6 ,e 6 ], [e 7 ,e 5 ], [e 8 ,e 5 ], [e 9 ,e 5 ], [ei ,e 5 ], 

[e 7 , e 6 , e 6 ], [e 6 , e 7 , e 7 ], [e 8 , e 6 ], [eg, e 6 ], [eio, e 6 ], [e 8 , e 7 , e 7 ], [e 7 , e 8 , e 8 ], 

[e 9 ,e 7 ], [eio,e 7 ], [e 9 ,e 8 ,e 8 ], [e 8 ,e 9 ,e 9 ], [ei ,e 8 ], [ei ,e 9 ,e 9 ], [eg, e w , e w ], 

[e 3 ,e 2 ,ei,e 2 ], [e 4 , e 3 , e 2 , e 3 ], [e 5 , e 4 , e 3 , e 4 ], [e 6 , e 5 , e 4 , e 5 ], [e 7 , e 6 , e 5 , e 6 ], 

[e 8 , e 7 , e 6 , e 7 ], [eg, e 8 , e 7 , e 8 ], [e w , e 9 , e 8 , e 9 ]) 

The first example presents the free Lie ring of rank two. The second and 
third examples are found in (Caranti, Mattarei, Newman & Scoppola, 1994) and 
are of great importance from the point of view of thin Lie algebras. The fourth 
example is due to Jiirgen Wisliceny and has a very nice structure, as we shall 
see later. The fifth one is not else than the presentation of the positive part of 
the classical Lie algebra A w factored out by the ideal generated by the torsions 
appearing over Z. 

The structure of L 4 is very nice as we mentioned before. Having a look at 
its canonical nilpotent presentation one has 

Theorem 6.1 The underlying abelian group of the lower central factors o/L 4 
is of the form 

L i+\ /L i+2 ^ f Cj/2 ® Coo © Coo © Coo ifi is even, 

1 \ Coo © Coo © Coo © Coo © Coo tfiis odd. 
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for < i < 67. 

By computation, the following structure theorem is true, for L 5 . 

Theorem 6.2 L 5 is nil-potent of class 10. Its lower central factors are of the 
form 

(9) LI/V+ 1 = C 00 <B---®C 00 , 

where the number if direct factors is 11 — i for 1 < i < 10. 

The structure of L\ is well known. (Hall, Jr., 1959) Theorem 11.2.2 provides 
a formula to compute the dimension for a lower central factor of L\. The result 
of the computation coincides with the known dimensions. 

The following table contains the computational informations concerning those 
examples. 



example 


class computed 


nr. of gens. 


torsion rank 


CPU time (ms) 


Li 


10 


226 





811100 


L 2 


14 


89 


70 


66200 


L 3 


12 


81 


63 


33150 


U 


7 


51 


22 


41333 


L 5 


6 


45 





254016 



Since all Lie rings presented by those relations are graded wc used the graded 
approach to compute the presentation for them. That is summarized in the 
following table. 



example 


class computed 


nr. of gens. 


torsion rank 


CPU time (ms) 


Li 


10 


226 





8633 


Li 


12 


757 





257383 


L 2 


14 


57 


38 


3150 


L 2 


18 


195 


171 


197483 


U 


12 


81 


63 


4133 


u 


67 


300 


32 


903553 


L 5 


10 


55 





36083 



All computations, except for the 67th class of L4, were done on a 486-DX2 
100 MHz PC with 16 MB memory plus 16 MB swap space. The above mentioned 
exception was computed on a Sun Sparccenter 2000 with 256 MB of RAM. 

One easily sees that using those ideas speeds up the computation with a 
significant factor. The need of using better algorithm for Smith normal form 
computation is reasonable. While computing L 3 an integer overflow occured 
with 32 bit arithmetic at the 68th factor. 
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